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Our aim is to understand how nature processes information. 
We can detect and describe the complex computational 
structure in natural processes and we can provide a 
quantitative characterization of essential aspects of sound 
generation in the flute which is the typical features of effect 
of turbulence inside the cavity. Aim of the study to establish 
the relation between the formations of particular sounds with 
inlet velocity and grove profile. This research is driven by 
the ever increasing amount of experimental data as 
delivered, to find internal effect cause to generate the sound 
together with change in the velocity and grove profile. 
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INTRODUCTION 

Nature has a very efficient way to process information. Processes studied by natural 
scientists involve systems that are either continuous, stochastic, spatially extended, or any 
combination of these, and fall strictly outside the range of discrete computation theory. The 
study of information processing in complex dynamical multiscale systems is, therefore, still 
in its infancy. This work is inspired from the wide complexity of the physical systems and 
consequently by the necessity to simplify their complexity into fundamental processes. 

In recent years, the Lattice Boltzmann Method (LBM) has developed into an 
alternative and promising numerical scheme for simulating fluid flows and modelling 
physics in fluids. The scheme is particularly successful in fluid flow applications involving 
interfacial dynamics and complex boundaries. Unlike conventional numerical schemes 
based on discretization of macroscopic continuum equations, the lattice Boltzmann method 
is based on microscopic models and mesoscopic kinetic equations. The fundamental idea of 
the LBM is to construct simplified kinetic models that incorporate the essential physics of 
microscopic or mesoscopic processes so that the macroscopic averaged properties obey the 
desired macroscopic equations. 
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Historically originating from the lattice gas automata (LGA) introduced by Frisch, 
Hasslacher, and Pomeau the lattice Boltzmann equation (LBE) has recently become an 
alternativemethod for computational fluid dynamics. The essential ingredients in any 
Lattice Boltzmann Models whichare required to be completely specified are: (i) a discrete 
lattice space on which fluid particles reside; (ii) aset of discrete velocities to represent 
particle advection; and (iii) a set of rules for the redistribution of particles residing on a 
node to simulate collision processes ina real fluid. Fluid-boundary interactions are usually 
approximated by simple reflections of the particles by solid interfaces. 

In a hydrodynamic simulation by using the lattice Boltzmann equation, one solves 
the evolution equationsof the distribution functions of pretended fluid particles colliding 
and moving synchronously on a highlysymmetric lattice space. The highly symmetric lattice 
space is a result of the discretization of particlevelocity space and the condition for 
synchronous motions That is, the discretizations of time and particle phase space are 
coherently coupled together. This makes the evolution of the lattice Boltzmann equation 
very simple, it consists in only two steps: collision and advection. One immediate limitation 
of the LBE method is due to its use of highly symmetric regular lattice mesh, which is 
usually triangular or square lattices in two dimensions and cubic in three dimensions. 
Obviously this is a serious obstacle to its applications in many areas of computational fluid 
dynamics. To deal with complex computational domains, various proposals have been made 
to use grids that are better suited to fit boundaries or to adapt meshes according to the 
physics of the system. 

It has been shown recently that the lattice Boltzmann equation is indeed a special 
finite difference form of the continuous Boltzmann equation with some drastic 
approximations tailored for hydrodynamic simulations ^ 5 > 7 b This makes the lattice 
Boltzmann method more amenable to incorporate body-fitted meshes 9] or grid refinement 
techniques [10 b In most cases the regular lattice mesh is abandoned by decoupling the 
spatial-temporal discretization and the discrete velocity set, so that interpolations can be 
used in addition to the advection on a non-regular or non-uniform mesh. 
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Figure II: Flute Size and their Classification 
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Figure III: Geometrical Representation of Cavity Generation 

In this work Lattice Boltzmann Equation is used to simulate typical flow of 
airthrough cavity. 

Flutes are unique woodwinds in their origin of sound. In order to produce sound, the 
pressure produced by the player’s mouth must be above atmospheric pressure. Thus, in 
order to create sound, a player must use power (work/time). When derived, this is 
equivalent to a mass times acceleration (Force) times a change in displacement over time 
(P= (m a Ad)/t). The acceleration of the air stream results in higher or lower pitches. The 
displacement refers to how far the air stream travels down the flute without leaving 
through any of its holes. In order to make music, a constant supply of energy is given by the 
player, because energy is lost to friction on the wood or metal of the flute, and some leaves 
through its holes. This causes resonance, which determines the frequencies of the notes. In 
a vertical flute, one side of the mouthpiece is sharper than the other. The air stream strikes 
against it, creating vibrations. Vibrations are made by a second hole below the mouthpiece. 
The bottom of this hole is filed to a 45° angle to create resonance. The velocity of air affects 
the sound created. For example, a quicker velocity created a high pitch. The velocity 
increases with an increase in pressure. Pressure equals force per given area (p = F/,4). The 

force increases as a person blows harder, and the area decreases as a person tightens 
his/her lips. Thus, the pressure is increased. 

To calculate the vibration of the air stream, one cycle is equal to the time for one 
pulse to travel over twice the total length of the flute. This accounts for both directions. The 
speed of sound in this case is “v,” so the frequency is equal to the speed of sound divided by 

twice the length [77 J. In terms of harmonics, their frequencies are in increasing order in 
terms of “n.” For example, the second harmonic has a frequency equal to ( 77 ), the third 

,... The picture above shows the first harmonic series on a flute tuned to C. Because its 

harmonics are in an even ratio (1:2:3:4:5:6:7:8:9:10), it sounds even and pleasing to the ear. 
This is called consonance. 
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Note 

Name 

M Length (Units 

of Octave) 

Factor to Multiply 
with Base Frequency 

Frequency 
Generated by Note 


C 

Sa 

1 

1 

65.4 

D 

Re 

9/8 

1.125 

72.56 

E 

Ga 

5/4 

1.25 

81.75 

F 

Ma 

4/3 

1.33 

86.982 

G 

Pa 

3/2 

1.5 

96.75 

A 

Dha 

5/3 

1.66 

108.56 

B 

Ni 

15/8 

1.875 

122.625 


Table 1: Octave Frequency Table 


THEORETICAL BACKGROUND 


In this Lattice Boltzmann model we use the kinetic theory as a fundamental of 
Lattice Boltzmann equation. In kinetic theory we consider the hard sphere model with 
elastic collision. The position and the momentum of the molecules are defined as x = (x,y, z) 

and p = (p.^p yJ p_J. Then f(x,p,t) is the probability density function for the presence of 

molecule in six — dimensional phase space. 

The probable no of molecules with position coordinates in the rangex x: + dxand 
momentum coordinates p 1 —> p + dp is given by f(x,p,t)dxdp. We introduce an external 


force F that is small relative intermolecular forces.If there are no collisions, then at 
time t -F dt, the new positions of molecules starting at x are x — (“T J dt = x + dx and the 

new momentum is p + Fdt = p — (~^) dt = P + dp. Thus, when the position and momentum 

are known at a particular time t, incrementing them allows us to determine at a future 
time t -f dt. Then the updated particle velocity distribution function is defined 

asf(x+ dx, p + dp,,t + dt). If the collision in not occurs then the particle velocity is not 
change then f(x -f dx, p -j- dp, t -f dt) — ffx,p, t) =0 and if the collision occur then there is a 

change in velocity distribution function and it is described 
as f(x T - dx, p + dp, t — dt) — f(x, p, t) =± 0. The Boltzmann Equation Expresses a balance 


between transport of a particle and collision between particles then 

Df Df 


Dt 


Tr-c nsonrt 


D t £ r[ igi 

df df 

Where C(f j models the pair wise collision between particles. 

C(f) is evaluated by the Tailor’s expansion of the updated velocity distribution 
function 

df i df m df (dxj 2 d 2 f (.dp} 2 d 2 f (at') 2 d 2 f 




+ ' 


2 dx 2 2 dp 2 2 
dxdp d z f dpdt d 2 f dtdx d 2 f 


2 dxdp 2 dpdt 2 dtdx 
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df ,0/ . df (dx) 2 d 2 f (.dp) 2 d 2 f (dt) 2 d 2 f 

dx—— 4- dv — 4- dt — 4--V 4- — : -V 4--— 

dx "dp at 2 dx 2 2 dp 2 2 dt 2 

dxdp d 2 f dpdt d 2 f dtdx d 2 f f — f*** 

2 dxdp 2 dpdt 2 dtdx r 

Then the equilibrium function is defined as f** = pw { |l -f -hc { u -f t^-C c,nj : - T^r^ 2 ] 

4 3 2 



Figure IV: The Lattice Pattern of the D2Q9 Model 

Wherew*, the weight factor corresponds to set of particles for D2Q9 lattice for rest 

4 1 

particle w 0 = 7 , nearest neighbours have w t — - for i = 1., 3,5,7 and farthest neighbours 
having w t = 77 for i = 2 ,4, 6 , 8 . c .is the discrete sets of velocities for D2Q9 lattice and is 
defined as c c = 0, nearest neighbours c ; = c | cos(i— l)- , sin(i — l) 7 ^ for i = 1 , 3 , 5/7 and 
farthest neighbours c t — \2c | cos(i — 1 )- , sinf; — 1)7 for i — 2,4, 6 , 8 . We take the speed 
of light in the media c = *J3w c 7 where Ax and At is the lattice spacing and time spacing 

respectively. For these models we take single time relaxation parameterr = 1 . 

The macroscopic variables are defined as a velocity moments of the velocity 
distribution function 

• Mass density: p(x, t) —m J f(x,p, t) d s p 

• Flow velocity: u(x , t) = — J p fix, p, t) d 2 p 


Temperature: T = 4- J Ip I ; f ix, p, t) d l p 


• Stress Tensor: p a ^ = f(x, p, t) (p a - u $ ) f(x, p,t)d s p 

SIMULATION RESULTS 


Note 


Velocity distribution 


Inlet /outlet average velocity 
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Table 2: simulation of different note of flute at time t= 4000 


Note 

Name 

Note 

Input 

Frequency 

Flow 

Velocity 

Output 

Frequency 

C 

Sa 

45.23 

75 

65.394 

D 

Re 

49.28 

82 

72.42 

E 

Ga 

55.55 

93 

81.85 

F 

Ma 

61 

101.5 

87.108 

G 

Pa 

65.20 

108.3 

96.613 

A 

Dha 

71 

118 

108.72 

B 

Ni 

83 

137.8 

121.14 


Table 3: Octave Frequency Table for Simulated Geometry 



■ Input Frequency 

■ Simulate 
Frequency 

■ Referance 
Frequency 
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CONCLUSION 

We simulate seven notes of flute as a diverse cavity models as depicted in figure III. 
Colum 2 in table 2 shows the turbulence generated by the different nodes and that origins 
for the particular sound in flute. Also in column 3 shows the average frequency at inlet 
(blue) and outlet (black) of the cavity. In table 3 demonstrate the inlet and outlet frequency 
generated by the notes that is a good argument with the simulated frequency and the 
reference frequency that can be shown in the above chart. 

In this work we consider equilibrium function in the Lattice Boltzmann equation is 
considered up to the second order with can be good arguments between the accuracy of the 
solution and computational needs. 
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